Predictions for high energy neutrino cross-sections from the ZEUS global PDF fits 
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We have updated predictions for high energy neutrino and antineutrino charged current cross- 
sections within the conventional DGLAP formalism of NLO QCD using a modern PDF fit to HERA 
data, which also accounts in a systematic way for PDF uncertainties deriving from both model 
uncertainties and from the experimental uncertainties of the input data sets. Furthermore the 
PDFs are determined using an improved treatment of heavy quark thresholds. A measurement of 
the neutrino cross-section much below these predictions would signal the need for extension of the 
conventional formalism as in BFKL resummation, or even gluon recombination effects as in the 
Ch . colour glass condensate model. 
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\ I. INTRODUCTION 

<N 

Predictions of neutrino cross-sections at high energies have sizeable uncertainties which derive largely from the 
' measurement uncertainties on the parton distribution functions (PDFs) of the nucleon. In the framework of the 
^p" 1 , quark-parton model, high energy scattering accesses very large values of Q 2 , the invariant mass of the exchanged 
O i' vector boson, and very small values of Bjorken x, the fraction of the momentum of the incoming nucleon taken by 
the struck quark. Thus when evaluating uncertainties on high energy neutrino cross-sections it is important to use 
the most up to date information from the experiments at HERA, which have accessed the lowest- a; and highest Q 2 
scales to date. The present paper uses the formalism of the ZEUS-S global PDF fits [l[ , updated to include all the 
C\J ■ HERA-I data. 

!> ' Conventional PDF fits use the Ncxt-to-lcading-order (NLO) Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
formalism 0, 0, 0, 0] °f QCD to make predictions for deep inelastic scattering (DIS) cross-sections of leptons on 
hadrons. At low- a; where the gluon density is rising rapidly it is probably necessary to go beyond the DGLAP 
formalism in order to sum ln(l/x) diagrams, as in the Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism [1, 0, Q 
(for recent work see 0, [ToL nTfl). An alternative approach is to consider non- linear terms which describe gluon 
recombination as in the colour glass condensate model [12| which has had considerable success in explaining RHIC 
data A recent suggestion is to use a structure function consistent with HERA data that saturates the Froissart 
unitarity bound and thus predicts a In 2 s dependence of the cross-section [l4| . Such approaches are beyond the scope 
of the present paper, which is concerned with the more modest goal of estimating the uncertainties on high energy 
neutrino cross-sections which are compatible with the conventional NLO DGLAP formalism. The motivation is to 
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provide an update on the neutrino cross-sections in the literature [lfj] which are widely used e.g. for estimating event 
rates in neutrino telescopes such as Baikal [l6j], ANTARES fl7l an d IceCub e [T el. cosmic ray observatories such as 
HiRes PJ] and Auger [20[, and radio detectors such as GLUE [2l|, FORTE [22^ RICE [H and ANITA [H- As a 
corollary, if cross-sections much outside these limits are observed, it would be a clear signal of the need for extensions 
to conventional formalism. To date no unambiguous signals which require such extensions have been observed. The 
prospect for measuring the cross-section using very high energy cosmic neutrinos in order to distinguish between 
theoretical suggestions for gluon dynamics at low x has been discussed by us elsewhere (2o| . 

Previous work on estimating high energy neutrino cross-sections (l5j used PDF sets which no longer fit modern 
data from HERA [26] and an ad hoc procedure for estimating PDF uncertainties. The present work improves on 
this in several respects. Firstly, we use a recent PDF analysis which includes data from all HERA-I running 
Secondly, we take a consistent approach to PDF uncertainties — both model uncertainties and, more importantly, 
the uncertainties which derive from the correlated systematic errors of the input data sets [27} • Thirdly, we use NLO 
rather than LO calculations throughout. Fourthly, we use a general- mass variable flavour number scheme p8l, [29j to 
treat heavy quark thresholds. 
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II. FORMALISM 

Partem Density Function (PDF) determinations are global fits [l|, HO, H2, HI, HH , which use inclusive cross-section 
data and structure function measurements from deep inelastic lepton hadron scattering (DIS) data. The kinematics 
of lepton hadron scattering is described in terms of the variables Q 2 , Bjorken x, and y which measures the energy 
transfer between the lepton and hadron systems. 

The double differential charged current (CC) cross-section for neutrino and antineutrino production on isoscalar 
nucleon targets are given by [35j |. 

dx dQ 2 ~ 2n(Q* + M^fx aMv)N) (1) 
where the reduced cross-sections a T {y(v)N) are given by 

a T {vN) = [Y + F»(x, Q 2 ) - y 2 F£(x, Q 2 ) + Y_xF£(x, Q 2 )] , (2) 

and 

a r (DN) = [Y + Fi(x, Q 2 ) - y 2 Fl(x, Q 2 ) - Y_xF*{x, Q 2 )] , (3) 

where the structure functions F2, xF^ and F^ are related directly to quark momentum distributions. 

The QCD predictions for these structure functions are obtained by solving the DGLAP evolution equations at NLO 
in the MS scheme with the renormalisation and factorization scales both chosen to be Q 2 . These equations yield the 
PDFs at all values of Q 2 provided these distributions have been input as functions of x at some input scale Qq. The 
resulting PDFs are then convoluted with coefficient functions, in order to obtain the structure functions. 

We use the PDF fit formalism of the published ZEUS-S global PDF analysis [l|, but this fit is updated as follows. 
First, the range of the calculation has been extended up to Q 2 = 10 12 GeV 2 and down to x = 10~ 12 . Second, all 
inclusive cross-section data for neutral and charged current reactions from ZEUS HERA-I running (1994-2000) are 
included in the fit. Third, the parametrization is extended from 11 to 13 free parameters, input at Q\ — 7 GeV 2 . 
In summary, the PDFs for u valence quarks (xu v (x)), d valence quarks (xd v (x)), total sea quarks (xS(x)), and the 
gluon (xg(x)), are each parametrized by the form 

Pl x P2 (l -x) P3 P(x), (4) 

where P(x) = 1 + P5X. The strong coupling constant is taken to be a s (M|) = 0.118 [36j. The total sea contribution 
is xS = 2x(u + d + s + c + 6), where q = q sca for each flavour, u — u v + M sca , d = d v + d sca and q — q sca for all other 
flavours. The flavour structure of the light quark sea allows for the violation of the Gottfried sum rule such that 
x{d — u) is non-zero, but only the normalisation of this quantity is free, the shape being fixed in accordance with E866 
Drell-Yan data [37| . A suppression of the strange sea with respect to the non-strange sea of a factor of 2 at Qq is also 
imposed, consistent with neutrino induced dimuon data from CCFR [38j |. The normalisation parameters, p\, for the 
d and u valence and for the gluon are constrained to impose the number sum-rules and momentum sum-rule. The 
low-x shape parameters P2 for the u and d valence quarks are set equal. Finally there are 13 free PDF parameters. 
Reasonable variations of these assumptions about the input parametrization are included in the model uncertainties 
on the output PDFs. 

A more important source of uncertainties on the PDFs comes from the experimental uncertainties on the input 
data. The PDFs are presented with full accounting for uncertainties from correlated systematic errors (as well as 
from statistical and uncorrelated sources) using the conservative OFFSET method. The uncertainty bands should be 
regarded as 68% confidence limits. A full discussion of approaches to estimating PDF uncertainties is given in [l| .l27| . 
The PDF uncertainties from this updated ZEUS-S-13 fit are com par able to those on the published ZEUS-S fit [l(, as 
well as the most recent fits of the CTEQ [H, Hi] and MRST [H 13] groups. 

Previous work [l5j treated heavy quark production by using a zero-mass variable flavour number scheme, with slow- 
rescaling at the b to t threshold. Although, as explained in Section [TTTJ the exact treatment of the b — > t threshold is 
not very important for the estimation of high energy neutrino cross-sections, it is important to use a correct treatment 
of heavy quark thresholds when determining the PDFs. We note that the central values of the sea quark distributions 
of the most recent CTEQ6.5 analysis [34[ which uses a general mass variable flavour scheme for heavy quarks, lie 
outside the 90% CL uncertainty estimates of the previous CTEQ6.1 analysis & which used a zero-mass variable 
flavour number scheme (as did all previous CTEQ analyses). This difference is significant for lower Q 2 (< 5000 GeV 2 ) 
and middling x (5 x 10 -5 <i<5x 10 -2 ) and this is a kinematic region of relevance to the present study. The heavy 
quark production scheme used in the present fit is the general mass variable flavour number scheme of Roberts and 
Thorne [H,[2gj]. The central values of both the CTEQ6.5 PDF analysis [H and the MRST2004 NLO analysis [3l| 
(which also uses a general-mass variable flavour number scheme) lie within, or very close to, the uncertainty bands of 
the present analysis over the entire kinematic region of interest. 
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FIG. 1: The PDFs and their fractional uncertainties at various Q 2 are shown for sea quarks (left) and gluons (right). 



III. RESULTS 

At very small x and high Q 2 , the vN cross-section is dominated by sea quarks produced by gluon splitting g — » qq. 
In this kinematic region, the parametrisation of the gluon momentum distribution is approximately: xg(x, Q 2 ) oc 
where A ~ 0.3 — 0.4. Figure [T] shows the predicted sea and gluon distributions from the present PDF fit and their 
fractional uncertainties, at various Q 2 values. This illustrates that the PDF uncertainties are largest at low Q 2 and 
at low-x. PDF uncertainties are also large at very high-x but this kinematic region is not important for scattering of 
high energy neutrinos. 

In QCD at leading order, the longitudinal structure function is identically zero, and the structure functions F 2 
and xF% for neutrino interactions on isoscalar targets can be identified with quark distributions as follows: 



F% = x(u + d+ 2s + 2b + u + d + 2c), xF£ = x(u + d + 2s + 2b - u - d - 2c). 
Similarly for antineutrino interactions, 



Fo 



x(u + d + 2c + u + d + 2s + 2b), xFg = x(u + d + 2c 



d-2s-2b). 



(5) 
(6) 



Assuming, s — s, c — c, b — b, we obtain F% — F% , whereas xF£ — xF£ = 2(s + s + b + b — c — c) = 4(s + b — c). At 
NLO these expressions must be convoluted with appropriate coefficient functions (such that F^ is no longer zero) but 
these expressions still give us a good idea of the dominant contributions. Note however that the b contribution will 
be very suppressed until Q 2 ^> M 2 ~ 3 x 10 4 GeV 2 , since the b — > t coupling is dominant. 
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FIG. 2: Predictions for F% , F£ and xF£ using a zero mass variable flavour number scheme, without the b contribution (upper 
panels), and with the b contribution (lower panels). 



In Figure [2] we show predictions for the neutrino structure functions F^, and xF£ and in Figure [3] we show 
the antincutrino structure function xF% . In order to illustrate the potential impact of the b contribution, these were 
calculated using the coefficient functions of the zero-mass variable flavour number scheme with (lower panels) and 
without (upper panels) the b contribution. Note that the input PDFs were still determined using the general mass 
variable flavour number scheme. 

The predictions for Fi and Fi, are somewhat suppressed without the b contribution, as is expected since the 
contribution of b to F2 is at most 20%. However, the effect on xF% is much more dramatic. This can be understood 
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FIG. 3: Predictions for xF£ for antineutrinos using a zero mass variable flavour number scheme, without the b contribution 
(left panel), and with the b contribution (right panel). 



by considering the LO expressions 

xFg = x(u v + d v + 2(s - c) + 2b) (7) 

and 

xF% = x(u v + d v + 2(c -s) -2b). (8) 

At low- a; the valence contributions are close to zero, while the strange and charm sea are of opposite sign and nearly 
equal, such that xF% is nearly all b quark and xFg ~ —xFg. 

Even though there are dramatic differences in predictions for xFa with and without the b contribution, this does not 
lead to significant differences in the vN and vN cross-sections because, at low-x, xF% < F2/5, and the y dependence 
suppresses the contribution of xF$ further. The b contribution to the reduced cross-section integrated over y is always 
less than ~ 25%. However, in practice it is even more suppressed for CC processes since b is important only at 
higher Q 2 (^> M t 2 ) because of the t threshold. Furthermore, in the total cross-section, the W propagator suppresses 
the contribution of the kinematic region Q 2 ^> Mfy, such that contributions from Q 2 ^> M 2 are suppressed very 
significantly. 

Figure [5] shows the predictions for the reduced neutrino cross-sections as a function of x for various Q 2 values above 
and below Mjy and M 2 . These illustrations have been made in terms of the reduced cross-section in order that 
one can see how the structure functions (hence the PDFs) contribute to the total cross-section. These cross-sections 
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FIG. 4: Neutrino-nucleon reduced cross-sections for various Q 2 at s = 3.6 x 10 7 GeV 2 i.e. E v = 1.9 x 10 7 GeV (left panel), 
and s = 10 10 GeV 2 i.e. E v = 5.3 x 10 9 GeV (right panel). 



have been calculated using the coefficient functions of the general-mass variable flavour scheme and are given for two 
representative values of the neutrino energy: s — 3.6 x 10 7 GeV 2 (=^> E v = 1.9 x 10 7 GeV) and s — 10 10 GeV 2 
(=>■ E v = 5.3 x 10 9 GeV). We do not show the antineutrino cross-sections separately because these are very close 
to the neutrino cross-sections at high energy. This is because the dominant structure function is F% = and 
although, xF^ ~ — xF% , the structure function xF% contributes with opposite sign in the neutrino and antineutrino 
cross-sections such that the net contribution of xFs is the same. 

The restriction on the lowest value of x probed for each Q 2 value is explained by the effect of the kinematic cut-off, 
y < 1: since, x = Q 2 /sy, we must have, x > Q 2 /s. This kinematic cut-off ensures that higher Q 2 values do not 
probe very low- a; until the neutrino energies are very high indeed. These figures illustrate which regions of x and Q 2 
contribute most strongly to the reduced cross-section for the different neutrino energies. The dominant contributions 
come from 50 < Q 2 < 10 4 GeV 2 (where the exact region moves up gradually with s); the contribution of higher Q 2 
(Q 2 > Myy) is suppressed by the W-propagator. For the lower energy E v = 1.9 x 10 7 GeV, the important range is 
10~ 6 < x < 1(T 3 , while for the higher energy E v = 5.3 x 10 9 GeV, this moves down to 1CT 8 < x < 1(T 4 . 

The PDF uncertainties are large at low- a; and low Q 2 , but since the dominant contributions to the cross-section do 
not come from very low Q 2 values, the PDF uncertainty on the total neutrino cross-section is quite small even at the 
highest energies considered here: s = 10 12 GeV 2 . 

The total neutrino cross-sections are now obtained by integrating the predicted double differential cross-section 
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d 2 a/dxdy with no cuts on either kinematic variable. 1 These cross-sections are tabulated in TableUat various values 
of s between 10 7 and 10 12 GeV 2 , together with their uncertainties due to the PDFs, including both model uncertainties 
and the experimental uncertainties of the input data sets. 2 The trend of the PDF uncertainties can be understood 
by noting that as one moves to higher and higher neutrino energies one also moves to lower and lower x where the 
PDF uncertainties are increasing. These energies are relevant to e.g. the Auger experiment where cosmic neutrinos 
can be detected both as quasi- horizontal deeply penetrating air showers and (specifically v T s) as Earth-skimming 
tau showers [2(J. We have shown elsewhere [25[ that the ratio of these two classes of events is a diagnostic of the 
vN cross-section, independently of the (rather uncertain) cosmic neutrino flux. The latter determines the absolute 
rates — e.g. assuming that extragalactic sources of the observed ultrahigh energy cosmic rays generate a neutrino 
flux saturating the "Waxman-Bahcall bound" [3^], it would take 10 years of running with a 3 x 10 4 km 2 array to tell 
whether the vN cross-section is suppressed significantly below the (unscreened) Standard Model prediction. Proposed 
satellite-borne detectors such as EUSO and OWL would scan even larger areas and achieve the necessary acceptance 
within a few years of running [40| . 







s [GeV J 


a{v) [pb] 


PDF uncertainty 


10 7 


1252 


±3% 


2 x 10 7 


1665 


±3% 


5 x 10 7 


2391 


±3.5% 


10 s 


3100 


±4% 


2 x 10 8 


4022 


±4.5% 


5 x 10 8 


5596 


±5.5% 


10 9 


7135 


±6% 


2 x 10 9 


9082 


±6% 


5 x 10 9 


12333 


±6.5% 


10 10 


15456 


±7% 


2 x 10 10 


19379 


±7% 


5 x 10 10 


25789 


±8% 


10 11 


31865 


±8% 


2 x 10 11 


39434 


±9% 


5 x 10 11 


51635 


±12% 


10 12 


63088 


±14% 



TABLE I: Neutrino-nucleon total CC cross-section, with the associated PDF uncertainty, at high energies. 



Figure O compares our CC cross-section to the widely used leading-order calculation of Gandhi et al [15[ which they 
fitted as: (cr^g/pb) = 5.53(E u /GeV) 0363 for 10 7 < (E v /GeV) < 10 12 . The present results show a less steep rise of 
the cross-section at high energies, reflecting the fact that more recent HERA cross-section data display a less dramatic 
rise at low- a; than the early data which was used to calculate the CTEQ4-DIS PDFs. A power-law description is no 
longer appropriate over the whole range 10 7 < {E v /GeV) < 10 12 , instead the relation 



In 



T NLO 
J CC 

pb 



ln(10 



98.8 



hi 



GeV J 



-0.0964 



(9) 



fits the calculated cross-section to within ~ 10% (P. Mertsch, private communication). 

Neutrino telescopes such as IceCube are optimised to probe lower energies of order a TeV 18]. In this case the 
high- a; region becomes important and the neutrino and antineutrino cross-sections are different because the valence 
contribution to xF% is now significant. In Table HT1 we give both neutrino and antineutrino cross-sections for values of 
s between 10 2 and 10 7 GeV 2 , together with estimates of their PDF uncertainties. The onset of the linear dependence 
of the cross-section on s for s < can be seen and in Figure [5] we compare the calculated <j/E v with some available 
recent experimental measurements [4l|, |42j , from the compendium on the Durham-HEPDATA database [43] . The 



1 Experiments may have specific cuts on e.g. y and we are happy to provide the differential cross-section for use in simulation programmes. 

2 Note that these are somewhat smaller than as shown in our previous work |25j | since we have now carefully evaluated the effect of heavy 
quark thresholds on the DGLAP evolution, which had been added on previously as a systematic uncertainty. 
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FIG. 5: The total CC cross-section at ultra high energies for neutrinos (left) and antineutrinos (right) along with the ±lcr 
uncertainties (shaded band), compared with the previous calculation by Gandhi et al [l5| ]. 



agreement is quite good, given in particular that our predictions are made for Q 2 > 1 GeV 2 (since perturbative QCD 
cannot be used at lower values); for s < 100 GeV 2 , there can be contributions of 0(10%) from even lower values of Q 2 
which are not accounted for here. However there are also corrections for nuclear shadowing which ought to suppress 
the cross-section by a comparable amount. A recent discussion of such effects can be found in ref . flij . 

The trend of the PDF uncertainties can be understood as follows: as one moves to lower neutrino energies one 
moves out of the very low- a; region such that PDF uncertainties decrease. These uncertainties are smallest at 10 -2 < 
x < 10~ x , corresponding to s ~ 10 5 . Moving to yet lower neutrino energies brings us into the high- a; region where 
PDF uncertainties are larger again. 







s [GeV 2 ] 


a{y) [pb] 


PDF uncertainty 


a(p) [pb] 


PDF uncertainty 


10 2 


0.334 


±3% 


0.151 


±4% 


2 x 10 2 


0.676 


±2.5% 


0.327 


±3.5% 


5 x 10 2 


1.69 


±2.5% 


0.864 


±3.5% 


10 3 


3.32 


±2% 


1.76 


±3% 


2 x 10 3 


6.47 


±2% 


3.55 


±2.5% 


5 x 10 3 


15.0 


±2% 


8.67 


±2.5% 


10 4 


27.6 


±2% 


16.6 


±2.5% 


2 x 10 4 


47.0 


±2% 


30.8 


±2% 


5 x 10 4 


89.4 


±2% 


64.8 


±2% 


10 5 


138 


±1.5% 


107 


±1.5% 


2 x 10 s 


204 


±2% 


171 


±2% 


5 x 10 5 


326 


±2% 


293 


±2% 


10 6 


454 


±2% 


423 


±2% 


2 x 10 6 


628 


±2.5% 


600 


±2.5% 


5 x 10 6 


937 


±2.5% 


915 


±2.5% 



TABLE II: Total CC cross-section for neutrinos and antineutrinos with their associated uncertainties at medium energies. 
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FIG. 6: The total CC cross-sections at medium energies for neutrinos and antineutrinos compared with with the previous 
calculation by Gandhi et al [l5| and with selected experimental data. 



IV. CONCLUSIONS 



We have calculated the charged current neutrino cross-section at NLO in the Standard Model using the best 
available D1S data along with a careful estimate of the associated uncertainties. As mentioned earlier, there are 
further uncertainties associated with QCD effects at very low x which are not addressed in the DGLAP formalism. 
When x is sufficiently small that a s ln(l/a;) ~ 1, it is necessary to resum these large logarithms using the BFKL 
formalism. Whereas such calculations at leading-log suggest an even steeper rise of the gluon structure function at 
low x (which would imply a higher v — N cross-section), this rise is not so dramatic at next-to- leading- log; for a 
recent application of NLL BFKL resummation to deep inelastic scattering see [Hj]. Moreover both the DGLAP and 
the BFKL formalisms neglect non-linear screening effects due to gluon recombination which may lead to saturation of 
the gluon structure function. This has been modelled in the colour dipole framework in which DIS at low x is viewed 
as the interaction of the qq dipole to which the gauge bosons fluctuate. An unified BFKL/DGLAP calculation (4|| 
supplemented by estimates of screening and nuclear shadowing effects, predicts a decrease of the v — N cross-section by 
20 — 100% at very high energies E v ~ 10 8 — 10 12 GeV [13] ■ An alternative approach uses the colour glass condensate 
formalism [l2j and predicts a similar suppression when a dipole model |48| which fits data from RHIC is used [49| . 
The predicted cross-section is even lower [49| if a different dipole model [5(j developed to fit the HERA data is used 
and the gluon distribution is assumed to decrease for x < 10~ 5 . Other possibilities for the behaviour of the high 
energy v — N cross-section have also been discussed [H], [5l| . 

Detectors for UHE cosmic neutrinos would be able to probe such new physics if they can establish deviations from 
the perturbative DGLAP prediction. Hence we recommend our calculated values for estimation of the baseline event 
rates in neutrino telescopes and for use in event generators such as ANIS [HJ. While the expected neutrino fluxes 
(e.g. from the sources of the observed high energy cosmic rays) are rather uncertain, experiments can in principle 
exploit the different dependence on the cross-section of the rate of Earth-skimming and quasi- horizontal events [25j . 
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